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ABSTRACT 

Accurate determination of spatial soil moisture distribution and monitoring its temporal 
variation have a significant impact on the outcomes of hydrologic, ecologic and climatic 
models. Development of a successful remote sensing instrument for soil moisture relies on 
the accurate knowledge of the soil dielectric constant (e t oii) to its moisture content. Two 
existing methods for measurement of dielectric constant of soil at low and high frequencies 
are, respectively, the time domain reflectometry and the reflection coefficient measurement 
using an open-ended coaxial probe. The major shortcoming of these methods is the lack of 
accurate determination of the imaginary part of In this paper a microstrip ring resonator 
is proposed for the accurate measurement of soil dielectric constant. In this technique the 
microstrip ring resonator is placed in contact with soil medium and the real and imaginary 
parts of e,oii are determined from the changes in the resonant frequency and the quality factor 
of the resonator respectively. The solution of the electromagnetic problem is obtained using 
a hybrid approach based on the method of moments solution of the quasi-static formulation 
in conjunction with experimental data obtained from reference dielectric samples. Also a 
simple inversion algorithm for e aoi i = e' + jej.' based on regression analysis is obtained. It is 
shown that the wide dynamic range of the measured quantities provides excellent accuracy 
in the dielectric constant measurement. A prototype microstrip ring resonator at L-band 
is designed and measurements of soil with different moisture contents are presented and 
compared with other approaches. 
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1 Introduction 

Successful modeling and understanding of global change are dependent upon an improved 
understanding of climatic and hydrologic processes. Soil moisture and its spatial and tempo- 
ral variations play an important role in both climatic and hydrologic models. Radar remote 
sensing of soil moisture has been an ongoing area of research over the past two decades. 
Characterization of dependence of dielectric constant of soil to its moisture content is an 
important step in this endeavor. In this paper the theoretical and experimental aspects of a 
field-portable device for accurate measurement of both the real and imaginary parts of soil 
dielectric constant are discussed. 

Two commonly used methods for field measurement of dielectric constant of soil are 
the time domain reflectometry (TDR) and the reflection coefficient measurement using an 
open-ended coaxial probe. In the TDR method the dielectric constant of the soil medium 
is determined from the delay time of the reflected wave from the open ends of a two-wire 
transmission line buried inside a soil medium [1]. The time delay can only provide the 
information about the real part of dielectric constant d T . The accuracy of the TDR method 
is limited by the resolution of the TDR system. In the coaxial probe method [2], the complex 
dielectric constant is estimated from the measured complex reflection coefficient of the open- 
ended coaxial line pressed against the unknown material. The accuracy of this technique 
is limited by the small dynamic range of reflection coefficient as a function of d T . Another 
problem with this technique when applied to soil medium is the small contact area of the 
coaxial tip which may be comparable to the soil particles. It should also be noted that soil 
dielectric measurement using coaxial probe is very sensitive to the applied pressure which 
degrades the accuracy of this measurement technique. 

In this paper a microstrip ring resonator as a device for the dielectric measurement of 
soil is proposed to overcome the aforementioned problems with the existing devices. Ap- 
plication of ring resonators in characterizing the dielectric properties of materials is not a 
new idea [3, 4, 5]. However, these techniques are specialized for dielectric measurements of 
thin materials in a stripline configuration and/or only characterization of real part of low 
dielectric materials is considered. Since both the real and imaginary parts of the dielectric 
constant of soil have a wide dynamic range with respect to the soil moisture content and 
the measurement configuration is different from those reported in the literature, a careful 
theoretical and experimental study of the proposed problem is required. The dielectric mea- 
surement technique is very similar to the standard cavity measurement technique [6] with 
a subtle difference. Basically, the real and imaginary parts of the dielectric constant are to 
be inferred from the shift in the resonant frequency and the change in the quality factor 
of the resonator after placing the resonator in contact with the soil medium. The subtlety 
stems from the fact that the microstrip resonator can be regarded as a partially-filled res- 
onant cavity and therefore the simple relationship between the dielectric constant and the 
measured changes in the resonant frequency and the quality factor can no longer be applied. 
Although the parti ally-filled nature of the microstrip resonators complicates the retrieval 
of dielectric constant from the measured quantities, it is a desirable feature in the soil di- 
electric constant measurements. In the standard cavity measurement, if the loss tangent of 
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the dielectric material exceeds 0.05, it becomes practically impossible to accurately measure 
the resonant frequency (/ 0 ) and the quality factor (Q) of the loaded cavity. However, for 
the partially-filled resonators the changes in /o and Q are substantially less affected by the 
loading dielectric material. It is shown that with appropriate microstrip resonator desigu, 
dielectric constant of soil with moisture content as high as 40% can easily be measured. 

In Section 2 the theoretical analysis of the problem is given. In this analysis which 
will be referred to as the forward problem a numerical solution for calculating the resonant 
frequency and the quality factor of the resonator in terms of e' and e" is provided. The 
numerical solution is constructed based on a quasi-static formulation. Since the practical 
microstrip resonators are constructed on finite substrate, the effect of finite substrate width 
is also investigated. In Section 3 the inverse problem is considered. There, a simple algorithm 
for retrieval of e' and e" from the measured changes in /o and Q based on regressional analysis 
is developed. Finally in Section 4 experimental results are presented. It is shown that the 
accuracy of the quasi-static solution degrades as the dielectric constant of the unknown 
material increases. This phenomenon has also been observed in stripline resonators [3]. In 
this section the theoretical results based on the quasi-static analysis are slightly modified 
to agree with experimental data obtained from reference dielectric samples. To verify the 
accuracy of this hybrid model, dielectric constant of a sandy soil with different moisture 
contents were measured using a prototype microstrip ring resonator and the results are 
compared with those measured by three other independent methods. 

2 Theory 

A microstrip ring resonator is a simple transmission line resonator whose geometry is shown 
in Fig. 1. The resonator is excited by a transmission line through a capacitive coupling 
and at certain frequencies, depending on the electrical length (perimeter) of the resonator, 
a standing wave pattern forms around the circular path of the resonator. The maximum 
voltage of the standing wave occurs at the exciting point. The resonant frequencies corre- 
spond to a condition where the perimeter of the ring is an integer multiple of the guided 
wavelength, that is, 


An — n — 1,2,3,- * * 

n 

where d is the diameter of the ring and \ g is a function of the microstrip parameters w s , h, 
e 3 and t 2 as shown in Fig. 2. 

At resonant frequencies there exists a voltage maximum at away from the excitation 
point. By placing a capacitively coupled transmission line at this voltage maximum point, 
the field in the resonator can be probed to detect the resonant frequencies. Basically the 
transmission coefficient S ? i of the two-port resonator is measured as a function of frequency 
and the resonant frequencies are identified as frequencies for which IS 21 I is maximized. It 
should be noted that the coupling capacitors tend to lower the resonant frequency and these 
must be as small as possible for accurate measurement of the resonant frequency. Spectral 
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measurement of |S 21 | can also reveal the quality factor of the resonator which is a measure 
of power loss in the resonator. The dissipated power in the resonator includes the dielectric 
loss, the conductor loss and the radiation loss. For a given microstrip resonator with known 
substrate dielectric constant e„ strip width u>„ and substrate height h we are seeking an 
algorithm that would enable us to compute the complex dielectric constant of the half-space 
medium in contact with the ring resonator from the spectral measurement of |S 2 i |. To estab- 
lish the relationship between the complex dielectric constant of the unknown medium and 
the measured guided wavelength and the Q of the resonator, the quasi-static approximation 
for microstrip transmission lines is used. Strictly speaking microstrip structures cannot sup- 
port TEM waves, however, at low frequency where the substrate height and strip width are 
small compared to the wavelength, it is expected that the quasi-static approximation pro- 
vides accurate results [7]. The quasi-static approximations found in the literature are usually 
concerned with lossless dielectrics. In what follows a brief derivation of Poisson’s equation 
for lossy dielectrics is given from which an integral equation for the charge distribution on 
the microstrip will be obtained. Starting from Maxwell’s equations for time harmonic fields 

V x H = jut'E + aE, ( 1 ) 

V-Z) = p c + p imp -D = ^E (2) 

where p c is the conduction charge density, p imp is the impressed charge density and e' and 
a are the permittivity and conductivity of the medium respectively. The conduction cur- 
rent J = <rE can be related to the conduction charge density p c through the continuity 
relationship V*J= — jwp c and therefore 

( 3 ) 

Combining (3) and (2) and defining the complex permittivity as e=e!-^-, it can easilv be 
shown that w 


— Pimp- ( 4 ) 

For TEM waves the electric and magnetic fields are irrotational and therefore the electric 
field can be obtained from a complex potential function (f> through 


E=-V<f> 


which together with (4) results 


V 2 <f> 


Pimp 


(5) 


Assuming the ground plane and the strip of the microstrip are at potentials zero and Vo, 
(5) can be solved subject to the boundary conditions. In this case p imp =p( x ) is the induced 
charge distribution on the strip. If G(x,y,x',y') represents the Green’s function of the 
problem, then the complex potential function can be obtained from 


# x >y ) ~ J_^ p(x')G(x,y,x',y')dx'. 
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Using the Fourier transform technique the Green’s function for the microstrip problem is 
found to be 


G(x — 



e ^ cos(a(x — x')) 
(e 4 coth(aA) + e 2 )a 


( 6 ) 


where t a is the dielectric constant of the substrate and e 2 is the dielectric constant of the 
upper half-space. Since the potential function is known over the metalic strip, the integral 
equation for the induced charge density can be derived by evaluating the potential over the 
strip and is given by 



x', h , h)dx ' . 


( 7 ) 


The integral equation can be solved numerically using the method of moments and point 
matching technique. Subdividing the strip into N sufficiently small cells, (7) can be cast 
into a matrix equation 


ZR=V 


where R is the vector of unknown charge density and V is the excitation vector whose entries 
are constant Vq potentials. To evaluate entries of Z, the Green’s function given by (6) must 
be computed for different values of observation (x m ) and source (x n ) points. The decay 
rate of the integrand of (6) (when y = h) as a function of or is rather poor (the asymptotic 
behavior of the integrand for large values of a is a ^ t< 1 +ta ^ )- To improve the convergence rate 
the order of integrations in (7) can be changed and the entries of the matrix are found to be 


2 r°° sin(q^)cos(a(x ro - x n )) 

7 r Jo a 2 [e, coth(ah) -f e 2 ] ° 


( 8 ) 


where Ax = It is noted that the impedance matrix is Toeplitz and symmetric, that is, 
Zmn = Z pq when m — n = p — q and Z mn = Z nm . Thus integral in (8) must be evaluated 
only N times to fully characterize the impedance matrix. Once the impedance matrix is 
characterized the charge density can be obtained from 

R=2~ l V . 


The quantities of interest in the microstrip resonator are the electrical length and the quality 
factor of the resonator which are directly related to the line capacitance and conductance. 
The line capacitance is proportional to the total charge in phase with the voltage and the 
line conductance is proportional to the conduction current. The total complex charge on the 
strip represented by Q = Q r + jQi can be obtained from 

Q = f_l p(x)dx~V< 

2 i j 

from which the line admittance can be calculated and is given by 

Y = y = 

V 0 Vq , ■ 
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It can be shown that over a narrow range of frequency the microstrip line behaves like a 

capacitor in parallel with a resistor. The line capacitance and conductance thus are given 
by 


C = ££ MZ7?] F/m 

« J 


( 9 ) 


G = s/m. ( 10 ) 

• i 

For wideband applications a more complicated equivalent circuit must be considered. As- 
suming that the substrate material and the upper half-space medium are nonmagnetic, the 
line inductance becomes only a function of geometrical feature of the line and can be ob- 
tained from magnetostatic analysis. An alternative approach is to find the line capacitance 
of the air-filled microstrip line (Co) and find the line inductance from 


L = 


1 

CoV? 


where V p =3xl0 8 m/s is the phase velocity of the air-filled fine [7J. To check the validity 
and accuracy of the numerical method, the line capacitance computed by the method of 
moments with e 2 = 1.0 and e. = 6.15 is compared with the line capacitance computed using 
the conformal mapping technique and is shown in Fig. 3. 

The quality factor of the resonator is defined as the ratio of the total average energy 
stored in the resonator to the power dissipated in a cycle. At resonance W e = W m = -CV 2 
thus the quality factor due to the dielectric loss can be obtained from [6] 4 


Qd = u 


We + Wr, 

P 


U 


\cv 2 

\GV 2 


uC 
G ' 


( 11 ) 


Another issue of practical importance is the finiteness of the substrate. The theoretical 
analysis given above is valid when the microstrip ground plane is of infinite extent. In practice 
the resonator substrate is finite and its effect on the line capacitance and conductance per unit 
length derived for microstrip with infinite ground plane must be investigated. Fig. 2 shows a 
microstrip with finite substrate and ground plane in contact with a half-space homogeneous 
medium. Our objective here is to characterize the smallest substrate width w g so that the 
difference between the line capacitance per unit length in this case and that of the infinite 
substrate is negligible. In this case a simple integral equation for the surface charge density 

on the strip and ground plane is obtained. The substrate dielectric is replaced with the 
induced polarization charge 


Pp — — 1) [6(r — r c )h • E — V • E) (12) 

where r c is a position vector that specifies the contour of the substrate (C) and n is an 
outward unit vector normal to the substrate boundary. After applying the Green’s theorem 
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the integral equation for the surface charge density and the complex potential function is 
found to be 


)<f>(x,y ) = -co(c, - 1 ) J c v') * ~ ~ dc + J_ w ' /2 p»( x ') G ( x i y; °) dx ' 

+ [ #/ Pg(x')G{x,y\x' ,-h)dx' (13) 

J -Wall 


where G(x,y\ x',y') is the Green’s function of the half-space dielectric problems and is given 
by 


G{x,y\x',y') = 


1 


2xe 0 


in yf{x — x') 2 + (y - y') 2 - “tt ln _ z ') 2 + {y + y') 2 • ( 14 ) 


Once the integral equation is solved numerically the line capacitance and conductance per 
unit length can be obtained as shown before. Fig. 4 shows the line capacitance of a finite 
substrate microstrip line with t, = 6.15, h = 0.245cm, and w s = 0.37 cm as a function 
of Wg/w s for 62 = 1. It is shown that when w s /w t > 5 the formulation based on infinite 
substrate is accurate for computation of the line capacitance. 


3 Inversion Algorithm 

In this section a retrieval algorithm for the complex dielectric constant of a half-space dielec- 
tric medium in contact with a microstrip ring resonator is sought. Following the standard 
cavity dielectric measurement approach the real and imaginary parts of the effective dielec- 
tric constant of the loaded resonator can be obtained from which the complex dielectric 
constant of the half-space medium must be inferred. Basically, first the resonant frequencies 
of the loaded (resonator in contact with the dielectric medium) and unloaded (resonator in 
free space) resonators are found from the measured frequency spectrum of |5' 21 1 . Then noting 
that at resonance the electrical length of the resonator is equal to the guided wavelength, it 
can easily be shown that: 

where f l and f u are, respectively, the resonant frequencies of the loaded and unloaded 
resonators. Since the ring resonator is a transmission line resonator, the relationship between 
the effective dielectric constant and the line parameters can be obtained from: 

= uJlc(\ - j ^)- 



Assuming the substrate is lossless (G u = 0) and noting that the line inductance is not a 
function of dielectric loading, it can be shown that 


C' = C 


u tjft 
E e// 


(16) 


7 



( 17 ) 


G‘ 

u l C* 



e 


< 7 / 


Comparing (17) to (11) and using (16), the imaginary part of the effective dielectric constant 
of the loaded resonator can be obtained from 


t 


» l _ fe/7 

e// Q d ' 


( 18 ) 


Using (15) and (18) the real part of the effective dielectric constant of the loaded resonator 
e 'fJ m terms of measured ^ and Q d and t u tf} can be obtained. Once and are 
obtained the line capacitance and conductance can be derived from (16) and (17). 

Besides the dielectric loss, other factors such as radiation and conduction losses deter- 
mine the quality factor of the resonator. If Q u represents the quality factor of the unloaded 
resonator and Q d is the quality factor due to the dielectric loss, the measured quality factor 
of the resonator (Q m ) can be obtained from [6] 


iz-'k+h ( w > 

Equation (19) can be used to measure Q d when the radiation losses of the loaded resonator 
and unloaded resonator are the same or the radiation loss is negligible compared to the 
conductor and substrate losses. This can be verified by measuring the quality factor of the 
resonator when loaded with different lossless dielectrics. However, if the radiation loss is 
significant and varies with dielectric loading, (19) can be modified to 

1 _ 1 1 1 

oi _ o: + o; + o; < 20 > 

where Q r the quality factor due to radiation loss and its dependency to e' must be determined 
empirically. As will be shown later, for most soil conditions Q d << Q u and therefore 
characterization of Q r with respect to e 2 is not necessary. 

The last step in the inversion algorithm is the computation of 4 and cr 2 from the mea- 
sured C and G . In the previous section a numerical forward model for computation of the 
line capacitance and conductance in terms of e 2 and cr 2 was developed. Computation of the 
forward model for various e' and <r 2 reveals that C l and G l are smooth functions of e' and 
<t 2 . Over the region of interest for e 2 and cr 2 these functions are one-to-one and onto which 
implies that the inverse functions exist. Since C(4, <r 2 ) and G( 4 , <r 2 ) are gentle functions of 
6*2 and <r 2 , the inverse functions are also gentle functions of C and G. To demonstrate this 
fact, and cr 2 versus the line capacitance and conductance of a microstrip with t a = 6.15, 
h = 0.245cm, and w, = 0.37cm are calculated numerically and are shown in Figs. 5 and 
6. It is shown that the line capacitance is strongly dependent on e' 2 and less influenced by 
<r 2 as expected. Also the line conductance is most sensitive to <r 2 and to a lesser extent 

to €2 Therefore the inverse functions can simply be approximated by their Taylor series 
expansions, that is 


C2 «a 0 + a x C + a 2 G + a 3 C G + a 4 C 2 + a 5 G 2 + a 6 CG 2 + a 7 C 2 G + a 6 C 3 + a 9 G 3 
d" a io C 2 G 2 + a u CG 3 -f a\ 2 C 3 G + ai 3 C 4 + aj4<j 4 


( 21 ) 
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( 22 ) 


(72 ~6 0 + hC + 62G + 63CG + + b 5 G 2 -(- b§CG 2 4~ bjC 2 G 4* b^C 2 4- b$G 2 

+bi 0 C 2 G 2 + buCGr 3 + buC 3 G 4 - b\3C* 4 " buG*. 

The coefficients of the polynomials can be obtained from the forward model using at 
least fifteen pairs of ( C , G). However, in order to find better estimate of the coefficients much 
more data points over the region of interest should be used in a least-mean-square estimator. 
Once the coefficients are determined the measured C and G can directly be used in (21) and 
(22) to find the unknowns e' 2 and 02. 

4 Experimental Results 

In this section the design considerations for the development of a prototype microstrip ring 
resonator and some experimental results are given. Prompted by the need for a high quality 
factor transmission line resonator, the ring resonator as opposed to other types of transmis- 
sion line resonators, such as a half- wavelength straight line resonator, was chosen. Another 
attractive feature of this resonator in measuring the dielectric constant of soil is its rela- 
tively large contact area with the soil medium. The coupling of energy to the resonator 
is accomplished using an unconventional method. Ring resonators are usually excited by 
a microstrip line through a capacitive air gap as shown in Fig. 7. However, for dielectric 
measurement applications it is better to excite the resonator with an open-ended coaxial 
line placed right underneath of the conducting strip as shown in Fig. 1. In this method 
the coupling capacitance of the resonator is not affected by the dielectric constant of the 
unknown material and therefore the quality factor due to the resonator coupling (Q u ) re- 
mains unchanged. Another important design parameter is the resonator filling factor. This 
parameter signifies the percentage of the stored electric energy in the substrate to the overall 
stored energy. For example to measure very lossy dielectrics a large filling factor is required 
to keep most of the field lines inside the substrate. The filling factor is directly proportional 
to the line capacitance of the unloaded microstrip (inversely proportional to the unloaded 
characteristic impedance). 

A microstrip ring resonator with a moderate filling factor that would allow accurate 
dielectric measurement of soil with volumetric moisture contents varying from 0% to 30% 
was designed using trial and error. It was found that a line with characteristic impedance of 
5017 would satisfy the filling factor requirements. For the prototype design, a Duroid 6006 
substrate with e s = 6.15 and h = 0.254cm was used. The strip width and the ring perimeter 
for the required unloaded characteristic impedance of 50fl and resonant frequency of 1.25 
GHz were calculated to be w a = 0.37cm and i = 11.35cm respectively. The dimension of the 
finite substrate was chosen so that the constraint w g /w s > 5 is satisfied for every point on the 
circular ring of the resonator. To examine the accuracy of the measurement technique and 
the aforementioned inversion algorithm, lossless reference dielectric samples were measured. 
An HP 8720 Network Analyzer with synthesized source was used to measure the quality 
factors and resonant frequencies reported in this paper. Six thick reference dielectric slabs 
(Stycast) with nominal dielectric constants of 2.2, 3, 6, 10, 13, and 16 were measured using 
the prototype resonator. It was found that the measurement accuracy degrades as the 
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dielectric constant of the medium increases. 

The experiments indicate that the quasi-static approximation overestimates the line 
capacitance when the dielectric constant of the half-space medium is increased. To investigate 
this phenomenon, two methods were pursued: (1) experimental, and (2) numerical. In the 
experimental investigation the effect of surface waves and the quasi-TEM nature of the 
approximate solution were examined. It is known that when the substrate height h is large 
the structure becomes capable of supporting surface waves in the substrate. Three thinner 
substrates were used to construct ring resonators (all with characteristic impedance of 50D) 
and the measurements with reference samples were repeated. Same amount of discrepancies 
wm obtained independent of substrate thickness. Next we examined the effect of quasi- 
TEM approximation in our formulation. For this purpose, a resonator was made from 
uroid 6010 substrate with e, = 10 and used to measure the permittivity of the Stycast 
sample with dielectric constant 10. In this case the medium surrounding the microstrip is 
homogeneous and the fundamental mode of propagation is TEM. However, a discrepancy in 
the measured dielectric on same order as before was observed. Having obtained inconclusive 
results from our experimental investigation, we resorted to numerical techniques. First a full- 
wave analysis for calculation of the effective propagation constant in microstrip lines based 
on the method of moments was used [8], It was found that for the frequencies and microstrip 
dimensions used m this investigation the quasi-static solution provided very accurate results. 

^ P o r n °^ e n m was ana ly zed completely using a three-dimensional FDTD numerical code 
19J. After 20,000 iteration corresponding to 7 hours CPU time on a Cray super-computer it 
was found that the resonant frequency obtained by the FDTD method had more discrepancy 
with the measured one than that obtained from the quasi-static method. One possible reason 

for the larger error is that the time domain methods may not be very suitable for resonant 
structures. 

As a fi Qal resort we decided to rectify the problem by including an empirical correction 
factor m the quasi-static solution. Let us denote the line capacitances derived from the quasi- 
static solution and those measured from the inversion algorithm by C q and C m respectively 
big. 8 shows the discrepancy in the line capacitance (AC = C q - C m ) as a function of the six 
measured dielectric constants. Since the measured discrepancy is a smootn function of e 2 and 

no discrepancy exist when e 2 = 1 the following expression for the corrected line capacitance 
is obtamed: 

C c = C q — k(t 2 — 1), (23) 

where k = 2.46 x 10" 12 is calculated using a least-mean-square linear estimator. This 
corrected line capacitance is then used in (21 ) and (22) to derive the coefficients a x • . • a 14 and 

i 14 ‘ These coefficients are listed in Table 1 for a microstrip with e, = 615 h = 0 254cm 
and w s = 0.37 cm. ’ ' ’ 

Next the behavior of the dielectric constant of a sandy soil as a function of moisture 
content is studied. The prototype ring resonator is used to measure soil samples with vol- 
umetric moisture contents ranging from 0% to 30%. Fig. 9 shows typical responses of the 

r !r nat ° r , m CO “ tact Wlth soil samples having moisture contents ranging from 0% to 
30 /o. 1 hree other independent techniques were also used to measure the soil samples. These 
include the open-ended coaxial probe, cavity resonator, and air-filled coaxial line. The dielec- 
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trie measurement using air-filled coaxial line is based on reflection coefficient measurement 
of a short-circuited coaxial line filled with soil samples. In this method the middle section 
of a long hollow coaxial line is used as the sample holder as shown in Fig. 10. Measuring 
the reflection coefficient over a wide bandwidth ( B ) and using the time domain capability 
of the network analyzer, the reflected signal after traveling twice through the sample can be 
measured. The hollow section of the coaxial line must be longer than the spatial resolution 
(3 x 10 8 /2 B) of the system. Fig. 11 shows the steps involved in measuring the reflection 
coefficient of the short circuit for a 50fl air-line of length l = 50cm filled with d = 20 cm of 
dry sand. The measured complex reflection coefficient is related to the dielectric constant of 
the sample through 

S n = -(1 - R i ) 2 e- i2kol e- j ' lko( ^ ; - l)d , (24) 

where R — is the reflection coefficient at the interface of air and soil. By inverting (24) 

the real and imaginary parts of the soil dielectric are obtained. The accuracy of this method 
decreases as the imaginary part of the dielectric constant increases unless the sample length 
is decreased. 

Figs. 12 and 13 show the measured real and imaginary parts of the dielectric constants of 
the sandy soil with seven different moisture contents using the four independent measurement 
techniques. The cavity method could only be used for the 0% and 5% moisture conditions 
and shows an excellent agreement with the ring resonator measurement. For the real part 
of the dielectric constant the agreement among all four techniques is very good. However, it 
is shown that the coaxial probe overestimates the imaginary part of the dielectric constant 
as the moisture content increases. 

5 Conclusions 

Theoretical, numerical, and experimental aspects of a microstrip ring resonator based di- 
electric measurement device are discussed in this paper. The real and imaginary parts of 
dielectric constant of a homogeneous medium in contact with the ring resonator is calculated 
from the shift in the resonant frequency and the change in the quality factor of the resonator. 
The solution to the forward problem is obtained from a numerical analysis of the quasi-static 
solution. For the inverse problem simple algebraic expressions based on regression analysis of 
the forward problem are also obtained. Measurements of reference dielectric samples showed 
that the accuracy of the quasi-static solution decreases as the dielectric constant of the half- 
space medium increases. To identify the source of error, full-wave numerical analysis and 
extensive experiments were conducted. Our efforts in identifying the source of errors were 
inconclusive. The errors in the quasi-static solution are corrected empirically from which 
the inversion expressions are obtained. The validity of the ring resonator measurement tech- 
nique is demonstrated by comparing the measured dielectric constants of a sandy soil with 
different moisture contents with those obtained from three other measurement methods. It 
is shown that the measurement technique is both accurate and versatile. The ring resonator, 
because of its large contact area, is very suitable for dielectric measurement of soil and can 
easily be constructed into a field portable device. 
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Figure 1: Top view and side view of a ring resonator used in this investigation. 



Figure 2: Geometry of a microstrip line in contact with a half-space homogeneous dielectric 
medium. 
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Figure 3: Comparison between method of mements (solid line) and conformed 
method (dash line) in calculation of the line capacitance of a microstrip line with 


mapping 
e,=6.15. 
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Figure 4: Line capacitance of a finite substrate microstrip line with c,=6.15 and tn a /ft=1.46 
as a function of w g /w s . 
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Figure 5: Contour map of the real part of dielectric constant as a function of line conductance 
and capacitance for a microstrip line with e,=6.15 and w,/h=\A§. 
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Figure 6: Contour map of conductivity as a function of line conductance and capacitance 
for a microstrip line with e,=6.15 and w t /h=lA6. 



Figure 7: Top view and side view of a conventional ring resonator. 
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Figure 9: Sand moisture spectral response. 



Figure 10: The air- filled coaxial line sample holder. 
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(a) Time Domain 



(c) Frequency Domain 



(b) Gated Response in Time Domain 



(d) Frequency Domain 



Figure 11: Steps in the measurement of the short circuit reflection coefficient: (a) time 
domain response of the total reflection coefficient, (b) time domain gated response of the 
short circuit, (c) magnitude and (d) phase of the short circuit reflecton coefficient. 
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Figure 13: Different methods for measuring ejf of a sandy soil with various moisture contents. 
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